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We report on the initial results from a measurement of the anti-neutrino flux and spectrum at a 
distance of about 800 m from the three reactors of the Palo Verde Nuclear Generating Station using 
a segmented gadolinium-loaded scintillation detector. We find that the anti-neutrino flux agrees 
with that predicted in the absence of oscillations excluding at 90% CL v Q — z? x oscillations with 
Am 2 > 1.12 x 10" 3 eV 2 for maximal mixing and sin 26 > 0.21 for large Am . 
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Nuclear reactors have been used as intense sources of 
P c in experiments searching for neutrino oscillations [0. 
These experiments usually detect v c by the process 
v c + p — > n + e + , where the cross-section- weighted en- 
ergy spectrum of i> c , peaking at about 4 MeV, can be 
deduced from the measured e + spectrum. Any v e flux 
deficit or distortions of the D c energy spectrum would 
indicate oscillations. The low energy of reactor v c al- 
lows these experiments to reach very small mass param- 
eters, albeit with modest mixing-angle sensitivity. Past 
experiments Q with detectors at 50-100 m from a re- 
actor have explored the mass-parameter range down to 
10 -2 eV 2 . The work described here, and a similar exper- 
iment elsewhere ||, are the first long baseline (~1 km) 
searches, designed to explore the parameter range down 
to 10 -3 eV 2 as suggested by the early Kamiokande atmo- 
spheric neutrino anomaly Q . Although later results from 
Super-Kamiokande B (appeared while this work was in 



progress) seem to disfavor the 



channel, a direct 



experimental exploration amply motivated this work. 

The Palo Verde neutrino oscillation experiment is lo- 
cated at the Palo Verde Nuclear Generating Station near 
Phoenix, Arizona. The total thermal power from three 
identical pressurized water reactors is 11.6 GW. Two of 
the reactors are 890 m from the detector, while the third 
is at 750 m. Our detector is placed in a shallow under- 
ground site (32 mwe overburden), thus eliminating the 
hadronic component of cosmic radiation and reducing the 
muon flux by a factor of ~ 5. The fiducial mass, seg- 
mented to reject the remaining background, consists of 
11.3 tons of 0.1% Gd-loaded liquid scintillator contained 
in a 6 x 11 array of 9 m-long acrylic cells, as shown in 
Fig. [|. Each cell is viewed by two 5-inch photomultiplicr 
tubes, one at each end. A v c is identified by space- and 
time-correlated e + and n signals. Positrons deposit their 
energies in the scintillator and annihilate, yielding two 
511 keV 7's, giving a triple coincidence. Neutrons ther- 
malize and are captured in Gd, giving a 7-ray shower of 
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FIG. 1. Schematic view of the Palo Verde neutrino detector. 



8 MeV total energy. 

The Gd loading of the scintillator has two advantages: 
it reduces the neutron capture time from 170 fj,s (on pro- 
tons) to 30 /is and provides a high energy 7 shower to 
tag the neutron capture, resulting in a substantial back- 
ground reduction. Both the positron and the neutron are 
triggered by a triple coincidence requiring at least one cell 
above a "high" threshold set at about 600 keV (positron 
ionization or neutron shower core), and two cells above a 
"low" threshold set at about 40 keV (Compton scatter- 
ing from annihilation photons or neutron shower tails). 
The triple coincidences are required to be within a 3 x 5 
matrix anywhere in the detector. 

The central detector is surrounded by a 1 m water 
shield to moderate background neutrons produced by 
muons outside the detector and to absorb 7's from the 
laboratory walls. Outside the water tanks are 32 large 
liquid scintillator counters and two end-caps to veto cos- 
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mic muons. The rate of cosmic muons is approximately 
2 kHz. The pattern of muons traveling through veto 
chambers and their timing relative to the central detec- 
tor hits are recorded for subsequent off-line analysis. The 
central detector is equipped with a system of tubes that 
allows the insertion of calibration sources in the small 
spaces between cells. In addition, a set of blue LEDs 
and optical fibers can produce flashes of light inside each 
of the cells. In order to reduce natural radioactivity, all 
building materials for the detector arc carefully selected, 
including the aggregate (marble) used in the concrete of 
the underground laboratory. 

The detector geometry, materials and electromagnetic 
interactions are simulated using the package GEANT. 
Hadronic interactions are described by FLUKA and the 
low energy neutron transport is simulated by GCALOR. 
The inclusive gamma spectrum from neutron capture on 
Gd is specially modeled according to measurements [||. 
Light quenching effects are also included H. 

Since the ultimate sensitivity of the experiment relies 
on a disappearance measurement, precise knowledge of 
the detector efficiency and of the expected v e flux from 
the reactors is essential. 

The efficiency calibration is based upon a primary mea- 
surement performed a few times per year with a cali- 
brated 22 Na e + source and an Am-Be neutron source. 
The 22 Na source is placed in the calibration pipes and 
mimics the effects of the positron from the D e interac- 
tion by providing annihilation radiation and a 1.275 McV 
photon which simulates the c + ionization in the scintilla- 
tor. The source is placed at 18 positions in the detector 
deemed to be representative of different conditions. The 
neutron detection efficiency is measured by scanning the 
detector with the Am-Be source where the 4.4 MeV 7 
associated with the neutron emission is tagged with a 
miniaturized Nal(Tl) counter. 

Other calibrations, used to measure the detector en- 
ergy response, are performed using the Compton edges 
from 137 Cs, 65 Zn and 228 Th sources. The same Th source 
is also used more frequently to track the scintillator trans- 
parency. Weekly runs of the fiber-optic and LED flasher 
systems are used, respectively, to monitor the gain and 
linearity of photomultipliers and the timing/position re- 
lationship along the cells. 

Since the energy deposition of the 511 keV 7's in one 
cell has a sharp falling spectrum (Compton scattering) 
it is vital to have the lowest possible "low" thresholds 
in the trigger and to understand the behavior of such 
thresholds with great accuracy. This second task is com- 
plicated by the fact that the trigger uses voltage ampli- 
tudes, while only charge from integrating ADCs is avail- 
able off-line. For this reason our detector simulation in- 
cludes a detailed description of the signal development in 
time. This code correctly describes the shape of pulses 
taking into account scintillator light yield, attenuation 
length and de-excitation time; photomultiplier rise- and 



fall-time and gain; and event position along a cell. The 
simulation of the detector response to the 22 Na source 
correctly describes the 40 keV (600 keV) threshold posi- 
tion to within 1.4 keV (2.6 keV), resulting in an uncer- 
tainty on the positron (neutron) efficiency of 4% (3%). 

The D e flux and spectrum from a fission reactor and 
the v c + p — * 11 + c + cross section are well known 
and are calculated by tracking the 235 U, 238 U, 239 Pu, and 
241 Pu fission rates in the three plant reactors, taking into 
account both power level and fuel age. The uncertainty 
in the P reaction rate is less than 3%. 

The data presented here was collected in periods of 
67.3 days in 1998 and 134.4 days in 1999. During the 
98 (99) data taking one of the far (near) reactors was off 
for 31.3 (23.4) days. While a detailed description of the 
data analysis will be reported elsewhere |l , here we will 
outline the principles of the analysis and the results. 

Neutrino candidates were selected by requiring an ap- 
propriate pattern of energy to be present in the detec- 
tor for the positron- and the neutron-like parts of the 
events. In addition the two sub-events are required to 
occur closer than about 1 m from each other. 

At our depth the background to v c events consists 
of two types of events: uncorrelated hits from cosmic- 
rays and natural radioactivity and correlated ones from 
cosmic-muon-induced neutrons. The first type can 
be measured by studying the time difference between 
positron- like and neutron- like parts of an event. By re- 
quiring that the time lapse between the two sub-events 
t CD be 5 [is < t cn < 200 [is, the uncorrelated background 
is reduced to 1.9 ± 0.1 events d _1 (3.7 ± 0.2 events d _1 ) 
for 1998 (1999), as measured from a fit to the capture 
exponential and a constant. 

The distribution of time intervals between a cosmic- 
ray (i crossing the detector and a z> c -like event carries 
information on the correlated background. From an ex- 
ponential fit we infer that the majority of correlated back- 
ground is produced by pairs of neutrons, where the cap- 
ture of the first neutron in each pair mimics the positron 
signature. The requirement that no cosmic-ray hits be 
present in a window of 150 [is preceding the 9 e candidate 
completes the event selection. While the detector effi- 
ciency, 77, is dependent upon the neutrino energy, the ef- 
ficiency integrated over the neutrino spectrum produced 
by the reactor (i.e. in the case of no oscillation) is ~ 8% 
(~ 11%) for 1998 (1999) (the higher efficiency in 99 is 
due to improvements in the data acquisition dead-time 
and trigger efficiency). 

The resulting rates N of v c candidates per day in 
different periods are given in Tab. |j| Along with the 
v e events this final data set contains the random back- 
ground mentioned above and a substantial amount of cor- 
related background. Two independent techniques were 
used to estimate and subtract the background. The 
most straightforward method ("Method 1") relies on the 
changes of the v c signal when different reactors are turned 
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off. For simplicity we first directly subtract the efficiency- 
corrected rates r = N/rj for "off" periods from the "on" 
periods. For f998 (f999) we compare the efficiency- 
corrected rates of observed interactions 95.0 ± 18.6 d _1 
(76.8 ± 13.6 d" 1 ) with the prediction of 63 d" 1 (88 d" 1 ). 
While the uncertainties quoted here are statistical only, 
we find good agreement with the hypothesis of no os- 
cillations. Wc can then proceed to add the 1998 and 
1999 data sets (which have a different combination of 
baselines) and obtain a neutrino energy spectrum for the 
difference "on" -"off". This spectrum well matches the 
model for no oscillations as shown in Fig. |^. 

In order to test quantitatively the oscillation hypothe- 
sis in the Am 2 — sin 2 29 plane we perform a x 2 -analysis 
comparing the expected to the measured rates. Informa- 
tion on the reactor neutrino flux, in the form of burn-up 
dependent fission rates, is computed for every detector 
run (on average ~12 h of data). This analysis effectively 
unfolds the background from the data using the i> -flux 
variation due to reactor power changes and fuel burn- 
up. At each point of the Am 2 — sin 2 28 grid the energy 
spectrum and the detector efficiency are calculated and 
used to predict the expected rates pi, where i denotes the 
different detector runs. We then construct the function: 
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The x 2 is minimized with respect to the background con- 
tribution b to the candidate neutrino rate Ni and the pa- 
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FIG. 2. Spectrum of energy deposited for neutrino events 
after "on" -"off" background subtraction. The sum of the 1998 
and 1999 distributions is compared with the Monte Carlo sim- 
ulation assuming no oscillations. Uncertainties are statistical 
only. The x 2 / n -d.f. between simulation and data is 8.5/8, 
supporting the no-oscillation scenario. 



Period 


98 "on" 


98 "off"* 


99 "on" 


99 "off"* 


Duration (d) 
V 

iV i^u ) 


36.0 
0.0746 
S8 9 + 1 


31.3 
0.0772 

"¥> 9 + 1 


111.0 
0.112 

52 9 + 7 


23.4 
0.111 
43.9 ± 1.4 


S v (d- 1 ) 
B (d" 1 ) 
i?Obs (d- 1 ) 

-RCalc (d _1 ) 


16.5 ± 1.4 
21.7 ± 1.0 
221 ± 18 
218 


13.4 ± 1.4 
18.8 ±1.0 
174 ± 17 
155 


25.2 ±0.9 
27.7 ±0.6 

225 ±8 
218 


15.1 ± 1.9 
28.8 ± 1.3 
136 ± 17 
130 



TABLE I. Summary of results from the Palo Verde exper- 
iment. The values in the second part of the table are derived 
from Method 2. B = B unc + B nn + B np . ifobs and 7?caic 
are the observed and calculated v e rates corrected by the ef- 
ficiencies r) for the case of no-oscillations. * Reactor at 890 m 
distance off. ' Reactor at 750 m distance off. Statistical un- 
certainties only. 



rameter a that accounts for possible global normalization 
effects. While statistical errors <Ji have to be individu- 
ally applied to each run, the systematic error <7 sys t is 
treated as a global parameter in the x 2 . The contribu- 
tion of different sources of systematics to cr sys t is given 
in the first column of Tab. ||. The effect of the z? e se- 
lection cuts is estimated by a random sampling of the 
unity hyper-volume defined by conservative ranges for 
each of the individual cuts. This technique properly takes 
into account the possible correlations between cuts. The 
90% CL acceptance region is defined in accordance with 
the procedure in (To) by Ax 2 > x 2 (Am 2 , sin 2 29) - Xb cst 
where x 2 (Am 2 , sin 2 26) describes the fit quality at the 
current grid point and Xbest that of the global best fit 
determined in the physically allowed parameter space. 
This procedure clearly prefers a no-oscillation scenario 
as illustrated in Fig. || (curve a). The fit also pro- 
vides an estimate of the background b = 19. 5 ± 1.3 d _1 
(6 = 26.3±1.7 d- 1 ) for 1998 (1999), and of a = 1.02±0.08 
with a x 2 /n.d.f. = 318/327. 

In an independent analysis ("Method 2") that will be 
described in detail elsewhere |llj we make use of the 
intrinsic symmetry of the dominant two-neutron back- 
ground to cancel most of the background directly from 
data and compute the remaining components from Monte 
Carlo simulations. This technique makes the best possi- 
ble use of the statistical power of all data collected. The 
rate of candidate events after all cuts can be written as 
N = B unc + B nn + B pn + S v where the contribution of 
the uncorrelated B unci two-neutron B nn and other corre- 
lated backgrounds B pn are explicitly represented, along 
with the v e signal S v . The dominant background B nn 
(along with B unc ) is symmetric under exchange of sub- 
events, so that an event selection with the requirements 
for the prompt and delayed event parts swapped will re- 
sult in a rate N' — B unc ± B nn + eiB pn + €2S V where 
ei and £2 account for the different efficiency for select- 
ing asymmetric events after the swap. We then calcu- 
late N — N' = (1 — ti)B pn + (1 - e 2 )S v where the effi- 
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Systematic 


Method 1 (%) 


Method 2 (%) 


e + efficiency 


4 


4 


n efficiency 


3 


3 


v c flux prediction 


3 


3 


v a selection cuts 


8 


4 


_B P n estimate 




4 


Total 


10 


8 



TABLE II. Origin and magnitude of systematic errors. Us- 
ing Method 2 for background subtraction reduces the system- 
atic uncertainty from the event selection cuts but introduces 
a new uncertainty due to the accuracy of the Monte Carlo 
used for the estimate of B pn . 



ciency correction £2 — 0.2 can be estimated from the D e 
Monte-Carlo simulation. We find that the processes of 
/i-spallation in the laboratory walls and capture of the 
/x's that are not tagged by the veto counter (4 ± 1)% 
contribute to (1 — ei)i? pn , while other backgrounds are 
negligible. Using Monte-Carlo simulation, we obtain 
(1 - ei)S pn = -0.9 ± 0.5 d" 1 (-1.3 ± 0.6 d" 1 ) for /Lt- 
spallation in 1998 (1999); the same figures for /i-capture 
are 0.6 ± 0.3 d" 1 (0.9 ± 0.5 d" 1 ) in 1998 (1999). This 
represents only a small correction to N — N' since ei is 
close to 1. While the Monte-Carlo model is accurate for 
the capture process, in the case of spallation we simulate 
the broad range of spectral indexes for the n-recoil energy 
reported in literature |i"T) . The average between differ- 
ent predictions is then used for B pn while the spread is 
used as an extra systematic error in the second column 
of Tab. [fj. Since no v c signal is present above 10 MeV, 
the observed integrated rate above such energy is used 
as a normalization of the Monte Carlo. The results are 
shown in the second part of Tab. [| for different running 
periods. Clearly Method 2 is also in agreement with the 
no-oscillation hypothesis. The excluded region, calcu- 
lated by comparing the expected and observed v e rates 
taking into account the effect of the oscillation parame- 
ters on r] and 62, is given in Fig. || (curve b). 

In conclusion, the data from the first period of run- 
ning from the Palo Verde detector shows no evidence 
for v c — P x oscillations. This result, together with the 
data already reported by Super-Kamiokande Q and 
Chooz II * excludes the channel — v c as being respon- 
sible for the atmospheric neutrino anomaly reported by 
Kamiokande Q . Data-taking at Palo Verde is scheduled 
to continue until the summer of 2000. 
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FIG. 3. Limits on mass difference and mixing angle from 
the present work (90% CL). Curves (a) and (b) are based on 
Method 1 and 2 for background subtraction, respectively, as 
described in the text. The method [10] is used to calculate 
the exclusion contours. The Kamiokande — v c atmospheric 
neutrino result is also shown. 
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